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ABSTRACT 

We present a new reduction pipeline for the VIRCAM@VISTA detector and describe the 
method developed to obtain high-precision astrometry with the VISTA Variables in the Via 
Lactea (VVV) data set. We derive an accurate geometric-distortion correction using as cali¬ 
bration field the globular cluster NGC 5139, and showed that we are able to reach a relative 
astrometric precision of about 8 mas per coordinate per exposure for well-measured stars over 
a field of view of more than 1 square degree. This geometric-distortion correction is made 
available to the community. As a test bed, we chose a field centered around the globular clus¬ 
ter NGC 6656 from the VVV archive and computed proper motions for the stars within. With 
45 epochs spread over four years, we show that we are able to achieve a precision of 1.4 
mas yr _1 and to isolate each population observed in the field (cluster. Bulge and Disk) using 
proper motions. We used proper-motion-selected field stars to measure the motion difference 
between Galactic disk and bulge stars. Our proper-motion measurements are consistent with 
UCAC4 and PPMXL, though our errors are much smaller. Models have still difficulties in 
reproducing the observations in this highly-reddened Galactic regions. 
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1 INTRODUCTION 

The VISTA Variables in the Via Lactea (VVV) variability cam¬ 
paign started in 2010. Thanks to the VIST A Infr aRed C amera (VIR- 
CAM, iDalton et al.l 120061 : lEmerson, McPherson. & Sutherland 
hood) , mounted at the 4.1 m telescope VISTA (Visible and Infrared 
Survey Telescope for Astronomy), this ongoing survey is mapping 
the Galactic bulge and disk to create a 3-D map of the Milky Way 
(Minniti et al. !20ld : ISaito et al.l2012l) . As for many long-term vari¬ 
ability surveys, the observing strategy is mainly focused on cover¬ 
ing a portion of the sky as large as possible in a single night, scan¬ 
ning the full field of view many times every few days. To this aim. 


* Based on observations with the 4 m VISTA ESO telescope, 
t E-mail: mattia.libralato@studenti.unipd.it 

% Visiting Ph.D. Student at STScI under the 2013 DDRF Spring program. 


the exposure time of each image has to be short enough in order 
to achieve the survey specifications. In the VVV survey, the typ¬ 
ical exposure time for A' s -filter images i s about 4 s, e .g ., a f actor 
7 smalle r tha n the 30-s threshold set by IPlatais et~aD 120021) and 
IPlatais. Wvse. & Zachariasl l l200d) as the minimum exposure time 
required to average out the large-scale semi-periodic and correlated 
atmospheric noise that harms ground-based astrometry. In spite of 
this, we chose to exploit the astrometric capabilities of this survey 
that will release to the community a data set with more than one 
hundred epochs over six years. 


In this paper, we present our reduction pipeline for the 
VIRCAM detectors and the geometric-distortion correction. As an 
example, we also show a few applications made possible by the 
astrometric accuracy reached by the VVV data set so far. 
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Table 1. List of the VIRCAM@VISTA data used for the astrometric cali¬ 
bration. Each observing block is made by A stC p images, where “step” is the 
dither spacing in arcmin between two consecutive exposures in an observ¬ 
ing block. The single-image exposure time is given by the integration time 
(DIT) multiplied by the total number of individual integrations (NDIT). 


Filter /V stcp Exposure Time 

Seeing 

Airmass 

(NDITxDIT) 

(arcsec) 

(sec z) 


Program ID: 488.L-05001A) - PI: Bellini 
NGC 5139 - w Cen 

February 23-24, 2012 

J 25 1.2 (6x10 s) 0.97-1.42 1.026-1.107 

J 258.5 (6x10 s) 0.74-1.08 1.134-1.198 


2 INSTRUMENT AND OBSERVATIONS 

VIRCAM is a mosaic of 4x4 detectors mounted at the focus of 
the VISTA 4.1 m telescope. Each detector is a Raytheon VIRGO 
2048x2048-pixel array and covers ~694x 694 arcsec 2 on th e sky . 
The average pixel scale is 0"339 pixeL 1 ^Sutherland et al.ll2014f) . 
The gaps between the detectors are quite large and correspond to 
42.5% and 90% of the detector size along the X and Y direction, 
respectively. 

Dithered observations are recommend ed to self-calibrate the 
geometric distortion of a detector (e.g., I Anderson et al J I2006L 
iBellini & Bedinl2010l fLi bralato e t al.l2014l) . However, the standard 
dither pattern adopted by VVV is not adequate for this purpose. For 
this reason, a calibration program (Program ID: 488.L-0500(A), PI: 
Bellini) was approved in 2012. The calibration field is centered on 
globular cluster NGC 5139 (cuCen). This field was chosen due to 
its relatively-high star density over more than one square degree. 

The field was observed in the J filter in two runs of 25 im¬ 
ages each (Table Q}, both organized in an array of 5x5 pointings, 
but with a different dither spacing (Fig. |T} Large dithers were 
taken to cover the gap between the 16 chips and to allow us to 
construct a single reference system for all observations. The small 
dither pattern was obtained to allow independent modeling of the 
high-frequency residuals of the geometric distortion for each chip. 


3 DATA REDUCTION 

We developed a reduction packa ge that makes use of the same 
tools described in lLibralato et al.l 1 20141) for the HAWK-I detector. 
Here, we briefly describe the software, and focus on the few 
differences between the two works. 

One peculiarity of VIRCAM is the striped pattern that affects 
all the images, both calibration and scientific. T hese stripe s are 
generated by the IRACE electronics ^Sutherland et alj[2014l) and 
change from one exposure to the next. To correct them, we made 
a FORTRAN routine based on the Cambridge Astronomy Survey 
Unit (CA SU) pipeline correctio n, that resembles the correction ap¬ 
plied bv lMavbhate et al.l d2008l) for the WFPC2@HST background 
streaks. We computed in each image the clipped-median value of 
the counts in each row, then we took the median of these values 
and subtracted it from the clipped-median value of each row. These 
differences represent the corrections to be applied to each row. We 



Figure 1. Outline of our adopted dither pattern used for 7-filter data. Large- 
and small-dither images are organized in two 5x5 arrays. The pointings of 
the large-dither images are taken to cover the gap between the 16 detectors. 
The total field of view covered by our observations is about l°7x2°0 on sky. 
The two zoomed-in panels show the dither spacing for large (azure panel) 
and small (red panel) dither pattern. 


did not include bad/warm/hot pixels while computing the median 
values. 

Using archival flat-field images we were not able to com¬ 
pletely correct for the pixel-to-pixel sensitivity variation of the de¬ 
tectors, in particular for chip [16] where the very high quantum- 
efficiency variation on short timescales sometimes makes it im¬ 
possible to properly apply the flat-field correctior[]. So, we con¬ 
structed master flat-field frames using the scientific images them¬ 
selves, masking all bad/warm/hot pixels and those in close prox¬ 
imity of any significant source (stars and galaxies) and considered 
these purged images as on-sky flat-field images. 

First, we applied dark and flat-field corrections to all im¬ 
ages. Then, for each chip we computed the median sky value in a 
5x5 grid and subtracted it according to the table. Then we made 
a 5x5 grid of fully-empirical PSF models for e ach dete ctor of 
each exposure, following the prescription given in I Anderson et al .1 
d2006h . Unlike the procedure given in the original img2psf_WFI 
and img2psf_HAWK-l programs, the finding criteria (minimum flux 
and minimum separation from brighter stars) to choose the stars 
that would be used to model the PSF are applied locally and are 
different in each cell of the grid. This way, we are able to find the 
most suitable combination of these criteria in each of the 5x5 re¬ 
gions of the chip (e.g., if the cluster center is located in a corner of a 
chip, the minimum separation from the brighter stars in that corner 
is usually lower than in the opposite corner where the crowding is 
lower). With an array of PSF models, we are able to measure po¬ 
sitions and fluxes for all sources on an image. The final catalogs 


1 http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/known-issues 
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(one for each chip) contain positions, instrumental magnitude^], 
and another quantity called quality-of-PSF-fit (QFIT) which repre- 
sents the absolute frac tional error in the PSF-model fit to the star 
dAnderson et alj|2008i) . The lower the QFIT, the better is the PSF 
fit. The QFIT parameter is a useful quantity to discriminate among 
well-measured and poorly-measured stars. Typically, in a> Cen cat¬ 
alogs we considered bright, unsaturated stars with a QFIT<0.05 to 
be well-measured stars. These selections allow us to always have at 
disposal over 100 stars per chip, with an average value of 350 well- 
measured stars for a corner chip and 1000 stars for a centermost 
chip. 


4 GEOMETRIC-DISTORTION CORRECTION 

In the large field of view (FoV) of VIRCAM, the tangential-plane 
projection effects are not negligible (at one degree from the tangent 
point this corresponds to more than 0.18 arcsec, ~0.5 VIRCAM 
pixel). This means the farther from the center, the larger is the dif¬ 
ference between the true position and the projected position of a 
star. 

We chose to perform an a uto-cal ibrat ion. By using as 
a reference sys tem 2MASS l lSkrutskie et al.1120061) or UCAC4 
dZacharias et alj|2013l) . which are among the most accurate ab¬ 
solute systems, we would have unavoidably ended up limited by 
their accuracy (of the order of 0.2-0.3 arcsec for 2MASS). Not 
to mention the non-negligible contribution from the stellar motion 
between the reference sy s tem a nd our exposures. Furthermore, as 
stated in lBellini & Bedir] d2010h . it is difficult to find a distortion- 
free reference frame with an homogeneous stellar density and lu¬ 
minosity. Therefore we adopted the auto-calibration solution. The 
basis of the auto-calibration is to observe the same star in as many 
different locations on the detector as possible and to compute its 
average position once it is transformed onto a common reference 
frame. Thanks to the large number and varied spacing of our dither 
pattern, a given star will be observed in several different locations 
in the FoV and, as such, the systematic errors in its mean posi¬ 
tion should average out. This way, the average positions of the stars 
should provide a reasonable approximation of their true positions 
in a distortion-free frame (the master frame). We built the master 
frame by cross-identifying the stars in each single-detector cata¬ 
log of each exposure. We used conformal transformations (four- 
parameter linear transformations: rigid shifts in the two coordi¬ 
nates, one rotation, and one change of scale) to bring the stellar 
positions measured in each image into the reference system of the 
master frame. In the left panels of Fig. [2] we show the effects of 
the projection on the master frame. The positional residuals along 
the X and Y axes show several bumps where two different chips 
overlap. 

When we first examined plots like these, it was clear that the 
bumps at the boundaries of the chips could be due either to internal 
distortions within each chip or to errors in placing the chips prop¬ 
erly with respect to each other. To ensure that the distortion within 
each chip was properly accounted for, we independe ntly solve for 
the ge ometric distortion of each chip, as done in iLibralato et al.l 
d2014l) . For this specific purpose, we only used the small-dither im¬ 
ages where the gaps between the chips are not covered (no detector 
overlaps with any other detector). The accuracy of the single-chip 

2 Defined as —2.5 X ]og(X counts), where 2 counts is the sum of the total 
counts under the fitted PSF. 


distortion solution was at the 0.02-pixel level (~7 mas). As an ex¬ 
ample, in the middle-left panels of Fig. [2] we show the residuals 
after we applied the single-chip correction to the catalogs and con¬ 
structed a small-dither-based master frame using only chip #10. As 
shown in Fig. [T| the small dithers do not allow us to put all the 
16 chips in the same reference system since the gaps are so large 
that the chips do not overlap each other. For this reason we selected 
one random chip (chip #10) to show that our single-chip distor¬ 
tion solution was good as we wrote above. Then, we applied our 
single-chip correction to all catalogs and used four-parameter lin¬ 
ear transformations to create a new master frame based on large- 
dither images. Again, the positional-residual bumps were still vis¬ 
ible (middle-right panels of Fig. [2}- Therefore, these trends in the 
positional residuals are ascribable to projection effects, which have 
to be taken into account while cross-identifying the catalogs. 

We chose to define a meta reference system in which to prop¬ 
erly project all single-chip catalogs and, at the same time, solve 
for most of the geometric distortion that affects this detector. We 
proceeded as follows. We used the 2MASS catalog as our initial 
reference frame. We projected the 2MASS catalog onto a tangent 
plane centered o n tuCen and followed the prescriptions given in 
Ivan de Ven et alj J2006I ) to convert R.A. and Dec. positions into 
pixel-based coordinates. This is an important step because we im¬ 
posed the master-frame scale to be exactly equal to 0.339 arcsec 
pi xel -1 for all chips . This value is the average pixel scale declared 
by [Sutherland.et ai] j20l4|) . 

Initially, we cross-identified all stars of each single-chip raw 
frames with the 2MASS catalog by using six-parameter linear 
transformations (which also include the deviation from orthogonal¬ 
ity and the change of relative scale between the two axes). Then, we 
located the center of each chip (x,y)=(1024,1024) on the 2MASS- 
based reference system. Without properly taking into account for 
the projection effects, the chip-center positions on the 2MASS ref¬ 
erence frame depend on their distances from the tangent point 
(at Cen center) and on the geometric distortion. To get rid of the first 
dependency and find the best position of the chip centers, we iter¬ 
atively de-projected the 2MASS catalog onto the celestial sphere, 
and then projected it again using as the tangent point the current 
chip-center position on the 2MASS reference system, in order to 
compute new, improved transformations. For each chip of each ex¬ 
posure/image we iterated the whole process five times (after the 
fifth iteration the adjustments were negligible). 

Once the chip centers in the 2MASS reference frame con¬ 
verged to fixed positions, in order to build the meta-reference sys¬ 
tem, we had to impose additional constraints. First, the meta center 
was defined as the average position of the four centermost chips. 
The second constraint we imposed is that the Y and X axes of our 
meta-frame system had to be oriented up and to the right, respec¬ 
tively. For each of the four centermost chips, we computed the an¬ 
gle between the expected meta-frame X axis and the segment that 
connects the center of the meta frame to the chip center. Then we 
rotated all chips by the average of the four angles. 

For each image we de-projected the 2MASS frame onto the 
celestial sphere and projected it back on a tangent plane, but this 
time using as tangent point the meta center computed as described 
above. Then, we rotated and shifted these 2MASS-based positions 
according to the other constraint. The final products of this ef¬ 
fort are 2MASS-based positions projected on the meta-frame cen¬ 
ter of each image, rotated and shifted to have the meta center in 
(jc,y)=(0,0). These positions represent the best approximation of the 
expected distortion-free meta positions. 

For each star in common between our catalogs and the 
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Figure 2. In each column of the Figure we show the positional residuals for different master frames. ( Left ): master frame based on single-chip, uncorrected 
catalogs and conformal transformations. ( Middle-left ): chip #10 master frame based on single-chip correction, small dithers, and conformal transformations. 
( Middle-right ): master frame based on single-chip correction, large dithers, and conformal transformations. Even if there is an improvement with respect to 
the “NO CORRECTION” case, the bumps (caused by projection-induced effects) are still visible. (Right): master frame based on our final distortion solution, 
six-parameter linear transformations and taking into account the projection effects. In the Top(Middle) row, we show the positional r.m.s. along the X(Y) axis 
as a function of the X(Y) position on the master frame. In the Middle-left and Middle-right-pane\ insets, we zoomed-in to better show the residuals. In the 
Bottom panels we plot the cr(Radial residuals) as a function of the J magnitude. The red dashed lines are set to the 3<x-clipped 68.27 th -percentile value of 
cr(Radial residuals). In all these panels we show only bright, unsaturated and well-measured stars. See text for more details. 


2MASS-based catalog, we have a pair of positional residuals that 
correspond to the difference between the raw-chip positions and 
the expected meta-frame positions (given by the stellar positions 
on the modified 2MASS reference frame). We used both saturated 
and unsaturated stars with magnitude 7<-12 and QFIT<0.2. Since 
2MASS is a shallow survey, we had to use saturated (7<-13.4) 


stars in order to have an adequate sample size. We divided each 
chip into a 3x3-grid elements and, in each such element, we 
defined the grid-point value as the average value of the residuals 
within. The cells have different sizes, with those close to the edges 
(for example 512x512 pixels on the corners) smaller than the cen¬ 
tral one (1024x1024 pixels), in order to better model the distortion 
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close to the edges. As described in lLibralato et alj l20l4) . for those 
cells adjoining the detector edges we shifted the grid points to the 
edge. We built a look-up table of correction for any location of 
the chip, using a bi-quadratic interpolation among the surrounding 
four grid points. To avoid extrapolation, our grid points extended 
to the corners, but this meant that we needed several iterations 
(each time applying 90% of the suggested correction to the raw 
positions and computing new residuals) before convergence could 
be achieved. 

After this first part of the correction, we have star positions 
transformed into a meta reference frame and corrected for geomet¬ 
ric distortion. All the stellar positions collected in one meta cata¬ 
log are those of the stars imaged in one exposure. Therefore, ac¬ 
cordingly to the Tabic |T| we have 50 meta catalogs at our disposal 
(25 of which are based on large-dither exposures, while the other 
25 are based on the small-dither exposures). The astrometric ac¬ 
curacy achieved is about 0.2-0.3 pixel, similar to that of 2MASS. 
The astrometric quality of our measurements should be ten times 
better than this, so to further improve our result, we applied an 
additional table-of-residuals correction to each chip by comparing 
the positions of the stars as measured in different meta catalogs, 
thus enabling a precision of ~0.03 pixel per comparison, as fol¬ 
lows. For each pair of meta catalogs (hereafter catalogs #1 and #2), 
we cross-identified all stars in common by using six-parameter lin¬ 
ear transformations. We found the meta center of catalog #1 into 
the reference system of catalog #2 and projected the stellar posi¬ 
tions measured in catalog #2 into the tangent plane centered at the 
center of catalog #1. Then we computed the positional residuals 
as the difference between the stellar positions in the meta #1 ref¬ 
erence system and the positions in the meta #2 reference system, 
once projected and transformed into the meta #1 reference system. 
For those meta catalogs obtained from the large-dither images, we 
compared each of them to the other 49 catalogs, while small-dither 
catalogs were only compared to the large-dither ones. When we 
computed the distortion correction for each chip individually, we 
used only small-dither images. We then applied this correction to 
large-dither images and looked at the residuals computed by com¬ 
paring our stellar positions with those of 2MASS. We noticed that 
the non-linear terms of the distortion over a very large scale were 
not completely accounted for. Therefore, we chose to compute the 
positional residuals by comparing only images far enough on the 
sky from each other. For each chip, we collected all these residuals 
together and divided them into an array of 1 lxl 1 square elements. 
We assigned to each array element the median value of the residu¬ 
als within. For any location on the chip, the correction is computed 
as the bi-linear interpolation between the surroundings four grid 
points. We iterated five times, computing new residuals and adding 
the new correction to the previous one. 

In summary, the distortion solution of each chip consists of 
two parts: a 3x3 look-up table of residuals (that is used to compute 
the correction at any inter-chip location via a bi-quadratic inter¬ 
polation between the surrounding four grid points), and an addi¬ 
tional fine-tuning 11x11 look-up table of residuals (this time using 
a bi-linear interpolation to compute the correction). The final stel¬ 
lar positions are distortion corrected and projected with respect to 
the center of the meta catalog. Therefore, each meta catalog is pro¬ 
jected into a different tangent plane. It is important to transform all 
the catalogs into the same tangent plane during the construction of 
the master frame. In the rightmost panels of Fig. [2] we show the re¬ 
sult of our efforts. We applied our distortion correction to the stars 
in each meta catalog. We used six-parameter linear transformations 


to bring these corrected positions on the master-frame reference 
system using the same tangent plane for each catalog. This way, 
the cr(Radial residualsf] improves from ~ 1.025 pixels (347.3 mas) 
to 0.023 pixel (7.9 mas). 

In Appendix[A]we show the distortion maps and the positional 
residuals along the X and Y axes, before and after the correction, 
for each of the 16 chips of VIRCAM. 

With this paper, we release a FORTRAN routine to correct 
the geometric distortion. It requires the single-chip raw coordi¬ 
nates (V aw ,y raw ) and the chip number. In output, the code com¬ 
putes (jc corr ,y corr ) coordinates in the meta-frame reference system. 
The code is available at our group’s web pagfj 


5 APPLICATION: NGC 6656 

To test our geometric-distortion correction we computed relative 
proper motions (PMs) of stars in the field of the globular clus- 
ter NGC 6656, (a, S) n <m.o = (18 h 36 m 23194. -23°54' I7'M. iHarrisI 
119961 2010 edition), using VVV data. We chose this object for its 
closeness and relatively-high PM with respect to the field objects. 
We used images taken between 2010 and 2014. We have 12 images 
in each of the 45 epochs used (except for one epoch for which we 
have 14 images) in the K s filter (from 2010 to 2014), while for the 
7 filter there are only 12 images taken in 2010. 

We obtained astro-photometric catalogs for each image of 
each epoch as described in Sect.[3j distortion corrected as described 
in Sect. [4] The VIRCAM photometry is calibrated by using stars 
in common with the 2MASS catalog. We applied linear relations 
between the VVV instrumental magnitudes and the 2MASS mag¬ 
nitudes based on well-measured, unsaturated stars. We replaced the 
photometry of VVV saturated stars with that of 2MASS. 

The adopted reference frame is based on images taken on Au¬ 
gust 16 th 2012 (which have the best available seeing, are the closest 
to the zenith, and are taken halfway between 2010 and 2014). The 
covered FoV is about l?lxl?5. We limited our PM analysis to the 
innermost region of the field, within a radius of 20 arcmin from 
the cluster center, where there is a significant number of cluster 
members. We then computed the coefficients of the local transfor¬ 
mations to tra nsform the stars’ posit ions of each image into the ref¬ 
erence frame dAnderson et alj|2006h . Local transformations reduce 
most of the uncorrected distortion residuals and other systematic 
effects that could harm our measurements. Indeed, the astrometric 
accuracy reached in our reference master frame is ~0.08 pixel (27 
mas), more than three times larger than that described in Sect. [4] 
The main reason for this larger uncertainty is that the VVV ob¬ 
servations are not taken with an astrometric strategy in mind (see 
discussion in Sect.[T]andO. Furthermore, our geometric-distortion 
correction is an average solution, suited for 7-filter images and at a 
specific epoch. The fact that the positional residuals are three times 

3 The cr(Radial residuals) is defined as: 


<x(Radial residual), = \j ——----——--- , 

T T . 

where {x. J . ,y f j) is the position of the i-th star of the j -th image, distor¬ 
tion corrected and transformed into the master-frame reference system, 
Tj is the the transformation of the j -th image into the master frame, and 
(^master, ymastery j s t ^ e distortion-free (master-frame) position of the i -th 
star. 

4 http://groups.dfa.unipd.it/ESPG/ 
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Figure 3. Multi-epoch fit of PMs for a field-like (“BULGE”) and a cluster¬ 
like (“M 22 HB”) star. (Top): in panels (a) and (c), the empty circles repre¬ 
sent star’s positions at different epochs transformed into the reference mas¬ 
ter frame, color-coded (as defined in the bottom panels) depending on the 
time interval relative to the reference epoch (J2012.62423). The master- 
frame position is represented by the solid black square (surrounded by an 
ellipse with semi-axes equal to the positional r.m.s. along the X and Y di¬ 
rection); the solid black circle is the expected position of the star at the 
reference epoch based on the PM fit. In the panel (a), the black arrow shows 
the ~4-yr displacement of the star. In panel (b) we show the position of the 
selected stars on the CMD. (Bottom): Motion in Y [panels (d) and (f)] and 
in X [panels (e) and (g)] as a function of the time from the reference epoch. 
The black line is the least-squares fit of the PM. 


larger also implies that the distortion correction is not stable over 
time scales of 6 months. The local-transformation approach com¬ 
pensates for these issues. 

In our local-transformation approach, we transformed the 
stellar positions as measured in each image into the reference- 
frame system using a subset of close-by, likely-cluster member 
(reference stars) to a given star to compute its linear-transformation 
coefficients. As such, our PMs are computed relative to the cluster 
mean motion, and cluster members will end up around the (0,0) 
location on the vector-point diagram (VPD). At the first iteration 
we selected the reference stars for the local transformation based 
on their position on the color-magnitude diagram (CMD). Once 
PMs were also estimated, we improved our reference-star list by 
removing all those stars which motion is not consistent with the 
cluster mean motion. 

We computed the stellar displacements as the difference be¬ 
tween the transformed single-exposure positions and the master- 
frame positions. In Fig. [3] we illustrate the multi-epoch PM fit for 
a Bulge star (left panels) and for a cluster member (right panels). 
In the top panels of the figure we show the stellar positions trans¬ 
formed into the reference frame, color-coded according to their 
epoch. In the bottom panels we show the displacements as a func¬ 
tion of time (relative to the reference epoch). The black lines are the 
weighted least-squares fit to the data, where the weight is defined 
as the square root of stars' QFIT (with poorly-measured stars’ hav¬ 
ing less weight). The proper motion along the two directions (ji x , 
p y ) is the slope of the straight lines. The relative PM errors are the 
formal errors of the least-squares fit. The constant terms xq and yo 
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Figure 4. (Left): Ks vs. (/— Ks) CMD of the NGC 6656 field. We show only 
well-measured-PM stars. We split the CMD in nine intervals of one mag 
each. The gray dashed line sets the average saturation threshold. The satura¬ 
tion level slightly varies within each VIRCAM chip. The variation becomes 
substantial across the total FoV, and of course from one exposure to the next 
(because of generally different seeing conditions). (Middle-left): VPDs for 
each of the corresponding magnitude interval. The mean motion of clus¬ 
ter members is centered at (0,0) in the VPDs. We plotted with red dots the 
cluster-like stars. The radius for the cluster-member selection (green circle) 
ranges from ~4.4 mas for stars with 16.5 <£s^ 17.5 to ~4.1 mas for those 
stars with 8.5<A!s<9.5. (Middle-right): Histograms for the /i a cos 8 proper 
motion distribution. The bin size changes depending on the total number of 
stars in each magnitude bin. Dual-Gaussian fit in black; individual Gaus- 
sians in red and azure are used for cluster and field n a cos 8 distributions 
respectively. The field distribution is wider than that of the cluster and con¬ 
taminate the cluster-member sample in all magnitude bins. (Right): CMD 
with only cluster-like-motion stars. It is clear that the fainter the magnitude 
bin, the higher is the field contamination in our sample. The PMs have been 
corrected for differential-chromatic refraction as described in the text. 


indicate the corrections to be applied to our reference-frame posi¬ 
tions at epoch J2012.62423. To exclude obvious outliers, for each 
star we iteratively removed one point at a time from the sample, fit 
the proper motion with the remaining points, and re-computed the 
vertical residual of the removed point. We rejected all those points 
for which vertical residuals were five times larger than the fit resid¬ 
ual. A final fit to all the remaining points provides our PM estimate 
for the star. 

In Fig. [4] we show the K$ vs. (J-K s ) CMDs and the VPDs 
for the stars in the NGC 6656 field. In the VPDs we show the 
(ji a cos 8,/is) PMs. In the left panel we plotted the CMD of the 
entire sample of well-measured-PM starqfl observed in the selected 
field, and split it into nine bins, one Ks -magnitude wide. In the 
middle-left panels we show the corresponding VPD for each mag- 


5 We plot the PM errors as a function of the Ks magnitude and drew by 
hand a fiducial line to remove obvious outliers. The cut is more impor¬ 
tant for faint stars, where PMs are less accurate. We used the same purging 
method for the stellar QFIT. Furthermore, we kept only those stars mea¬ 
sured in at least 50 exposures. 
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Figure 5. (Bottom-left): Ks vs. (J-Ks) CMD of NGC 6656. We plotted the 
horizontal-branch stars with light blue open squares and red-giant-branch 
stars with light red open circles. We considered only stars with Ks between 
13.5 and 15. (Top-left): VPD with only the blue and red stars before the 
differential-chromatic-refraction correction. The azure solid square repre¬ 
sents the mean motion of the blue stars, while the red solid circle that of 
the red stars. We also plot the 3 cr error bar as reference. (Bottom-right): 
p a cos 6 and pg as function of the (J-Ks) color. A linear fit (solid black 
line) in each direction was used to model the correction. (Top-right): VPD 
after the correction. 


nitude bin, in which we drew a circle to enclose cluster-like-PM 
stars. The radius is a compromise between including field stars 
with a cluster-like mean motion and excluding cluster members 
with larger PM uncertainties. Saturated and faint stars have 
poorly-measured PMs, so the selection radii are more generous. 
In the middle-right panels we show the histograms of the motion 
along the p a cos S direction (where the difference between cluster 
and field stars in the PMs is more evident) to show the field 
contamination expected by our selections. The rightmost panel 
shows the CMD of NGC 6656 members. The horizontal branch, 
red-giant branch, main sequence turn-off and upper main sequence 
are mostly cleaned by field stars and can be used for studies of 
the properties of these stars in the near infrared. We also analyzed 
the possible impact of differential-chromatic-refraction effects on 
our PMs (Fig. [5). We selected two samples of NGC 6656 stars 
with a Ks magnitude between 13.5 and 15, one on the horizontal 
branch (blue stars) and the other on the red giant branch (red stars). 
Even if the mean motion of both samples is the same within the 
error bars, we chose to remove the small contribution of this effect 
(clearly visible in the /r a cos<5 direction). All PMs in Fig. [4][6] and 
[7]have been corrected accordingly. 

To evaluate our PM precision we proceeded as follow. By 
construction, our PMs are relative to the cluster mean motion. 
As such, the mean location of cluster members on the VPD is 
(0,0), and the observed dispersion of cluster members should in 
principle reflect the stellar internal-motion dispersion plus our 
measurement errors. The internal dispersion of NGC 6656 is about 
0.5 mas yr~' (assuming a d istance of T2 kpc and a central velocity 
dispersion of 7.8 km s~ lHarrislfl996l. 2010 edition). The 1-D 


Average PM Bulge - Disk 

(A/^cosd.A/^j) = (-2.13±0.22,-3.58±0.24) mas yr -1 -> (4.18±0.32) mas yr -1 



Figure 6. (Top-left): Ks vs. (J-Ks) CMD. We plotted Disk main sequence 
stars with blue open circles and Bulge red-giant-branch stars with red open 
squares. We limited our samples to 12.8</fs^l4 in order to use only well- 
measured bright stars. We did not plot all stars within 5 mas yr~* from the 
center of the VPD to exclude most of the NGC 6656 members. (Bottom- 
right): VPD with stars within 12.8<AT^14. The green circle used to ex¬ 
clude cluster members has a radius of 5 mas yr -1 . (Bottom-left): histograms 
of the iug for the Bulge and Disk stars previously selected. We fitted each 
histogram with a single Gaussian. (Top-right): as on Bottom-left but for the 
p a cos 6. 


dispersion (defined as the 68.27 th percentile of the distribution 
around the median) of bright, unsaturated (13<X S <14) cluster stars 
in the VPD is of about 1.5 mas yr -1 . By subtracting in quadrature 
the internal dispersion of 0.5 mas yr -1 , we end up with an ex¬ 
ternal estimate of our PM precision, which is of about 1.4 mas yr~*. 

To further test of our astrometric accuracy we measured the 
relative difference between the Bulge and the Disk bulk motion 
within the same selected VVV field of NGC 6656 (Fig. |6). To 
this aim, we selected two samples of stars, one from the Disk 
main sequence and one from the Bulge red giant branch. We 
considered only Disk (Bulge) stars that in the Ks vs. (J-Ks) 
CMD are bluer (redder) than the respective fiducial line of the 
sequence. Furthermore, we considered only those stars with PMs 
larger than 5 mas yr -1 with respect to the bulk motion of the clus¬ 
ter. We fit a single Gaussian to the histograms of the Bulge and 
Disk PMs along each direction, and found a relative displacement 
of (A p a cos d,AyUj)=(-2.13±0.22,-3.58±0.24) mas yr -1 . The abso¬ 
lute difference between Bulge and Disk bulk motions is therefore 
4.18±0.32 mas yr~'. 

To test this result, we measured the relative displacement 
between the Bulge and the Disk components using the motion of 
the same test stars a s measured in the UCAC4 and the PPMXL 
feoeser, Demleitner, & Schilbachi 1201 (it) catalogs. We found a 
relative displacement of 3.79±0.98 mas yr~* using UCAC4 and 
2.93±1.3 mas yr~‘ using PPMXL, which are in agreement with our 
estimate within the error bars, though our estimate has a smaller 
uncertainty. We also compared our measured difference between 
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Bulge and Dis k mo tion with that predicted by the Besan§on 
models jRobin et alj|2003i) . We simulated both populations in the 
same field covered by our application. We adopted an exponential 
trend for photometric and PM errors as function of the magnitude 
to create a model as close as possible to our data. The major 
challenge was to take into account for the correct absorption 
toward the Galactic pla ne. We used the Bulge Ex tinction And 
Metallicity (BEAM, see iGonzalez et alj|20l3 . 1201 3f ) calculator to 
compute the average extinction in our field. This value, divided by 
NGC 6656 distance, gives us the diffuse absorption of 0.15 mag 
kpc~'. The difference between Bulge and Disk motion obtained 
this way is 1.38±0.12 mas yr -1 . This value is not consistent with 
our measurements. We performed different simulations varying 
the absorption coefficients to understand if the absorption law 
could somehow change the simulated kinematics, but we found the 
results were about the same. We attribute this significant difference 
to the difficulty of the Besanjon model in simulating the reddening. 
Galaxy stellar densities and kinematics toward the Galactic Plane 
where the extinction is high. 


5.1 Future perspectives 

The VISTA Variables in the Via Lactea will be completed in 2016 
dHempel et alj|2014l) . and the time baseline provided by the uni¬ 
form VVV data will be about six years. As an example, we com¬ 
bined the VVV images o f NGC 6656 and the HAWK-I data pre¬ 
viously used bv lLibralato et alj d2014t) . Since the HAWK-I images 
were taken in 2007, we used only the VVV archival images be¬ 
tween 2010 and 2013 in order to have approximately the 6 years 
of time baseline. We computed the PMs as described in the previ¬ 
ous section and in Fig. [7] we show the resulting CMDs and VPDs. 
As expected, with a larger time baseline we are able to completely 
separate cluster and field stars. This example shows again the great 
astrometric potential of the full-baseline VVV data. Older epochs 
(both optical and near-infrared data) are available in the archives, 
and the proper motions will be an invaluable resource to distinguish 
the different stellar populations in the Galaxy. 


6 CONCLUSIONS 

In this paper we present our reduction pipeline for the VIRCAM 
detector and the geometric distortion solution based on the J fil¬ 
ter. Thanks to our distortion correction and to the adopted dithered 
observing strategy, we are able to reach a positional residual of ~8 
mas in each coordinate in each exposure across the entire FoV of 
VIRCAM. Note that we are talking about relative astrometry. Our 
absolute astrometry is not as good as the relative one because the 
linear terms are constrained only with 2MASS. 

We release a FORTRAN routine to correct the geometric distor¬ 
tion. For a given position in a single-chip raw frame (V aw ,/ aw ) and 
the chip number, the code produces (v corr ,y corr ) coordinates in the 
meta-frame reference system. The code is available at our group’s 
web pag^H The use of this distortion solution is encouraged regard¬ 
less of the specific method adopted to measure stellar positions. 
Each meta catalog is projected into a plane tangential to its center. 
This offers the best single-catalog, distortion-free positions. Please 
note that, in order to construct a common reference frame, all meta 

6 http://groups.dfa.unipd.it/ESPG/ 
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Figure 7. (Top): vector-point diagrams with a time baseline of about six 
years obtained combining the HAWK-I and VVV data. (Bottom): Ks vs. 
(J-Ks) color-magnitude diagrams of the stars in common between the 
HAWK-I and VVV fields. We plot only well-measured-PM stars. On the 
left panels we plot the entire sample, while in the center and right panels 
we plot only cluster members and field stars respectively. We considered 
NGC 6656 stars those stars with a proper motion within 4 mas yr 1 around 
the cluster mean motion (red circle centered in (0,0) in the middle VPD); 
while the stars enclosed in the ellipse centered in (—11.5,1.2) mas yr 1 with 
major and minor axes of 12 and 9 mas yr~* are probable field stars. 


catalogs should be instead projected into the same tangent plane 
(see Sect.0. 

As a test bed of the astrometric accuracy reached by our 
geometric-distortion correction, we applied our reduction pipeline 
to a set of VVV archival images. We chose a field centered on the 
globular cluster NGC 6656 and we computed the relative proper 
motion of the NGC 6656 and Galactic bulge and disk stars, as well 
as the individual motion of each star in the field. We noticed that 
our astrometric accuracy is worse (~0.08 pixel) using VVV data. 
Our geometric-distortion correction is an average solution and the 
distortion is not entirely stable. However, by starting with a good 
average solution, local transformations (used to compute the proper 
motions) can be used to efficiently achieve optimal precision even 
with this type of data. We demonstrate that we are able to separate 
cluster and background/foreground field stars with a time baseline 
of only four years. The cluster stars, in the cleaned CMD, can be 
used for the study of the stellar populations of NGC 6656. We also 
showed that the field stars, in the direction of NGC 6656, are of 
great use, e.g., to separate (and study) the proper motions of the 
Galactic disk and bulge components. We demonstrated that our re¬ 
sults are consistent with what can be obtained using UCAC4 and 
PPMXL catalogs, though our measurements have a much smaller 
error. Galactic models fail to reproduce the observations, likely be- 
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cause of the difficulties to reproduce the reddening and kinematics 
towards the Galactic bulge. 

With the images analyzed in this paper and a time baseline 
of about four years, we obtained a typical astrometric precision of 
1.4 mas yr _1 for bright, unsaturated well-measured stars. This value 
corresponds to cr,.— 21 km s _1 at the distance of NGC 6656 (3.2 kpc 
from lHarrii 1 1 9961 . 2010 edition), or ~53 km s -1 at 8 kpc (a ref¬ 
erence distance for the Bulge). At the end of the VVV survey, the 
total time baseline will be of about six years, thus further increasing 
the final achievable PM accuracy. The use of older, archive, optical 
and near-infrared data will further enhance the proper-motion capa¬ 
bility of the VVV survey. The astrometric capability of this survey 
is complementary to GAIA, in particular in the most crowded and 
heavily-absorbed regions not reachable by GAIA, and to study ob¬ 
jects below its magnitude limit (G~20). 
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APPENDIX A: GEOMETRIC DISTORTION MAPS 

In this appendix we show the distortion maps (Fig. ED and the 
positional residuals (Fie, IA2l and Fig. [A3} for the 16 chips of VIR- 
CAM before and after we applied the distortion solution described 
in Sect. [4] 
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Figure Al. ( Left ): Residual trends for the 16 chips when we use uncorrected stellar positions. The labels on the top-right comer of each box represent the chip 
number. The size of the residual vectors is magnified by a factor of 250. Some degree of distortion is clearly visible in the outermost chips. (Right): Residuals 
after our distortion correction is applied. The size of the residual vectors is now magnified by a factor 5000. 
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Figure A2. From the Bottom-left panels, clockwise: 8x vs. X, 8y vs. Y, 6x vs. Y and Sy vs. X for each of the 16 VIRCAM chips before we applied the distortion 
correction. 
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Figure A3. As in Fig. lA2l but after the distortion correction is applied. 






































































































































































































































































































































































































